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Stability  Analysis  of  a  Rigid  Body  with  a  Flexible  Attachment 
Using  the  Energy-Casimir  Method. 
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Abstract 

We  consider  a  system  consisting  of  a  rigid  body  to  which  a  linear  extensible  shear  beam  is 
attached.  For  such  a  system  the  Energy* Casinur  method  can  be  used  to  investigate  the  stability  of 
the  equilibria.  In  the  case  we  consider,  it  can  be  shown  that  a  test  for  (formal)  stability  reduces 

to  checking  the  positive  definiteness  of  two  matrices  which  depend  on  the  parameters  of  the  system 

"  ^  » 

and  the  particular  equilibrium  about  which  the  stability  is  to  be  ascertained.  . 
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1.  Introduction 


We  consider  a  rigid  body  to  which  a  long,  flexible  appendage  is  attached.  A  coordinate 
reference  frame  is  fixed  in  the  rigid  body  with  the  origin  at  the  center  of  mass  of  the  rigid 
body.  The  flexible  attachment  is  assumed  to  lie  along  the  second  coordinate  asds  when  the 
configuration  is  at  rest.  (  see  figure  1.)  The  equations  of  motion  for  such  a  configuration,  under 
suitable  assumptions  and  with  the  appendage  modeled  as  a  linear  extensible  shear  beam,  are 
derived  by  Krishnaprasad  and  Marsden  in  [2].  In  deriving  the  equations  of  motion  they  use 
Hamiltonian  methods  in  the  context  of  Poisson  manifolds  and  reduction,  (see  [2]  for  the 
explicit  formula  for  the  Poisson  brackets  involved.) 

If  we  assume  that  the  momentum  of  the  system  which  arises  from  the  appendage  rotating 
with  the  rigid  body  is  negligible,  then  our  Hamiltonian  is  of  the  form 


H  =  J  J-p  .  p  +  { /' 11=^  +  i /'xg  .  g 


We  assume  that  J  is  the  inertia  matrix  of  the  rigid  body  and  that  po  is  uniform  mass  per 
unit  length  of  the  attached  appendage  of  length  t.  The  reduced  phase  space  b  coordinatized 
at  any  time  by  w,  the  convected  angular  velocity  vector  of  the  rigid  body;  r(s),  the  convected 
displacement  of  the  shear  beam  at  a  point  s,  0  <  s  <  ^;  and  ni(s)  the  momentum  density  of 
shear  beam  at  the  point  s.  The  vector  p  b  the  body  angular  momentum  vector  of  the  rigid 
body,  thus  p  =  Jw.  Finally,  K  b  the  diagonal  matrix  of  elastic  coefficients. 


In  our  investigation  we  are  interested  in  the  stability  of  the  system  about  equilibria  points. 
These  equilibria  will  satisfy, 


0  =  Ju)Xb>  +  ax  K— 
da 

^  1 

0  =  — m  +  r  X  w 
Po 

d^T 

0  =  +  m  x<j 

ds‘ 


(1.2) 

(1.3) 

(1.4) 


Two  boundary  values  are  associated  with  these  equations, 


In  [2],  a  stability  algorithm  based  on  the  Energy*  Casimlr  method  was  applied  to  a  spe¬ 
cific  family  of  equilibria  (see  section  4.2  below).  The  essence  of  the  stability  algorithm  is  to 
recognize  that  the  relevant  Poisson  structure  { « ,  *  }  admits  nontrivial  Casimirs  i.e.  functions 
F  that  Poi$$on~commute  with  any  function  of  the  phase  space.  It  follows  that  these  are  also 
conserved  quantities  for  the  dynamics  of  (1*1).  Specific  Casimirs  may  be  found  such  that 
the  relative  equilibria  defined  by  (1.2)  -  (1.4)  are  critical  points  of  (H  -I-  C^)  on  the  reduced 
phase  space.  Formal  atability  follows  from  establishing  definiteness  conditions  for  the  second 
variation  X?^(H  -t-  C^)  at  the  relative  equilibria.  To  establish  rigorous  nonlinear  atability,  one 
has  to  carry  out  certain  convexity  estimates  as  in  [2]. 

The  purpose  of  this  paper  is  to  establish  a  ayatematie  procedure  for  carrying  out  the  formal 
atability  atep  for  arbitrary  equilibria  satisfying  the  equations  (1.2)  -  (1.4).  This  has  useful 
applications  in  the  engineering  context  where  the  model  at  hand  represents  the  mechanics  of  a 
spinning  spacecraft  with  a  flexible  attachment  (such  as  a  boom  for  carrying  instruments  or  an 
antenna).  See  [2]  for  related  remarks  and  references.  The  procedure  derived  here  recovers  the 
results  of  [2]  when  applied  to  the  specific  example  considered  there,  (see  section  4.2  below.) 


2.  Computation  of  the  First  and  Second  Variations 

In  this  section  we  compute  the  first  and  second  variations  of  the  Hamiltonian  plus  the 
Casimir  function,  H  +  C^.  From  the  previous  definitions  of  these  we  know 

H  =  J  J-P .  p  +  J  /^‘  W  +  1  *  i..  (2.1) 

and  the  Casimir  function  may  be  taken  to  be 


f  rxmdsll*). 

Jo 


We  will  denote  the  first  and  second  variations  by  i>(H  +  C^),  and  +  C^).  Note  that 

because  of  the  distributed  nature  of  the  system  we  are  dealing  with  we  will  need  to  compute 
variational  derivatives  instead  of  ordinary  gradients. 

2.1.  Computation  of  the  First  Variation 

For  the  integrals  in  the  Hamiltonian  we  consider  variational  differentials  defined  by 


Thus,  letting 


Similarly,  let 


c-»0  €  Jq  Sx 

(2.3) 

Jo  Po 

(2.<) 

/•*  1 

Df2{m)=  1  — m'^mds. 

Jo  Po 

(2.5) 

(2.6) 

«r*- 

(2.7) 

If  we  integrate  this  equation  by  parts  with  the  boundary  conditions  ST(t)  =  ^r(O)  =  0,  we  get 

^  •  Sr  ds.  (2.8) 


For  the  integral  term  in  the  Casimir  function  we  are  taking  variational  derivatives  of  a  cross 
product  term.  If  we  define 


lrf  =  l|P+ /‘rxindsll*. 

Jo 

then 

Plla||’ =  2a  •  (tfp  +  f  Tx6md6+  f  Srxmds). 

Jo  Jo 

If  we  combine  all  of  the  above  we  get  the  expression  for  the  first  variation 

i?(H  +  C^)  =  J~*p  •  fp  +  /  — m*  finds  —  /  K-^-z  •  Srda 

Jo  fio  Jo 

+  ^'(||a||*)  a  •  (fp  +  f  Tx6mds+  f  Srxmds). 
Jo  Jo 


(2.9) 

(2.10) 


(2.11) 


2.2.  Computation  of  the  Second  Variation 

For  the  second  variation,  the  starting  point  is  the  expression  for  the  first  variation.  The 
terms  arising  from  the  original  Hamiltonian  are  straight  forward  to  compute,  they  are 


I?(J-»p.fp)  = 

J  *fp.fp. 

(2.12) 

1 

Dll  — m.fmds)  = 
Jo  Po 

f*  1 

/  — Sm-Smds, 

Jo  Po 

(2.13) 

f*  d^r 

(2.14) 

Note  that  we  can  use  the  boundary  conditions  on  Sr  to  get 

/o  ds  ds 

(2.15) 

Next  we  consider  the  component  which  arises  from  the  Casimir  function  which  we  added  to 
the  Hamiltonian.  IVom  the  first  factor  of  this  term  we  compute, 

I?^'(a)  =  2^"(l|a|l’)  a*  (fp+ /  rxSmds+f  Srxmds).  (2.16) 

Jo  Jo 

From  the  second  factor  of  the  Casimir  term  we  compute 

D{pi’{Sp+  f  rxSmds+f  Srxmds))  = 

Jo  Jo 

fjSp  +  I  T  X  Smds+  f  Srxm ds|j* 

Jo  Jo 

+  2(p  +  J  r  X  mds)  -  {J  Sr  x  Sm ds).  (2.17) 


We  use  the  above  to  get  the  expression  for  the  second  variation 


P*(H  +  Ci)=J  *5p*5p+  f  — 6m- 6m  da  +  f 

'  Jo  Po  Jo  da  da 

+  2^"(l|all*)  (a*  (5p+  f  t  x  6mda+  f  frxmds))* 

Jo  Jo 

+  ^'(l|a||*){||fp  +  J  T  X  6mda  +  j  6t  X  mdsll^ 

+  2(p  +  I>  X  m da)  -  (J  6t  X  6m  ds) }. 


(2.18) 


3.  Computation  of  a  Stability  Criterion 

The  conditions  which  assure  that  the  first  variation  2?(H  +  C^)  at  an  equilibrium  is  zero 


^'(lla*I|»)a*xr*  =  --m% 

Po 


^'(Kina*xm‘  =  -K:^, 


where  u*  =  J“^p*,  and 


a*  =  p*  -t-  /  r*  X  m*  da. 

Jq 


We  use  the  superscript  e  to  denote  evaluation  at  an  equilibrium.  If  we  dot  (3.1)  with  a*  we 


If  we  evaluate  the  first  variation  at  an  equilibrium,  incorporating  the  above,  then  we  can  derive 
conditions  which  assure  the  stability  of  the  equilibrium.  In  the  following  sequence  of  steps  we 
demonstrate  how  this  is  done. 

Step  1  :  Evaluate  the  Second  Variation  at  an  Equilibrium 

Recall  the  second  variation.  If  we  use  the  above  to  substitute  for  ^'(||tt*||^)  in  this  ex¬ 
pression  and  rearrange  slightly  we  find  that 

P’(H-f  C^)(p.^,Bi.)  =  J"'fp-fip-l-  j  ^6m-6mda  +  J 


w*-a* 


Step  S:  Collect  Terms  Containing  6p 

Now,  collect  together  terms  in  which  the  quantity  appears.  Our  expression  for  the 
second  variation  at  an  equilibrium  can  then  be  written 

i?*(H  +  C^)  =  •  fp  -  +  25p  •{[  r*  X  fmds  +  /  £t  x  xn‘ ds)) 

L  lla*ll  Jo  Jo 

-f- 2^"(||a*||*){(a*  •  fp)*  +  2(a*  •  5p)(a*  •  (^  t‘ x  Brads +  j  ^rxm*ds))}j 
-  — — ^11  f  T*  X  Bra  d8+  f  Bt  xra*  dsll* 

iia‘ir'yo  Jo  " 

+  2^"(l|a‘||’)(a*  {f  r‘  X  Binds  +  f  BTXm‘  ds))* 

Jo  Jo 

<l»*  •  tt*  €  /  r  r  J  ^  C  £  J  V  J 

-2 - ‘11  Br  X  Bmds)-\-  /  — Bra  -  Brads  +  /  K-^ — -^ds 

l|a*l|^  Vo  Jo  Po  Jo 

(3.10) 


Step  4‘  Complete  the  Square 

The  term  in  square  brackets  which  contains  the  Bp  terms  can  be  rewritten 
.  ]  =  (J-‘  -  +  2^ (||«*||>K  ®  «‘)<P  •  Sp 

+  2(- +  20"(]|a*|!*)a*  ®  ili*)fip  •  (  /  T‘xSmdB+f  Srxm‘ds) 

||a*ll  Jo  Jo 

(3.11) 

In  this  expression  we  use  ®  to  denote  the  tensor  product  and  I  the  identity.  Note  that  a*  0a* 
is  a  tensor  of  rank  2.  We  can  complete  the  square  for  this  expression  provided  the  quantity 

has  an  inverse. 

We  next  assume  this  inverse  exists  and  define  the  two  symmetric  matrices  M  and  N  by, 
M’-M  =  J-'  -  +  2,S"(||a'f  )a*  8  a' 

=  J7'  (SIS) 

N^'M  =  -  ^^1  +  2#"(1|«'||')<.'  8  a* 


i»,  i»- 


4»>.  *>,  Ata  •»■.■>«  *»- 


Completing  the  square  for  the  term  in  brackets  we  now  get 

•  =||M^p  +  N(/  T*xSxnds+f  5rxm*<fs)||* 

Jo  Jo 

—  N^N(  f  r*  X  Sm  d8+  f  Sr  x  m*  da)  ■  {  f  t*  x  Sm  d8+  f  St  xm*  da] 
Jo  Jo  Jo  Jo 

(3.15) 


The  term  in  braces  is  bounded  below  by  a  perfect  square  when  N^N  >  0.  For  this  to  be  the 
case  we  need  to  assume  that  the  inverted  matrix,  J7*  is  positive  definite,  in  general  it  need 
not  be.  Note  that  this  assumption  will  impose  conditions  on  ^"(||a‘|p).  The  requirements 
on  the  parameters  in  this  matrix  to  assure  it  is  strictly  positive  definite  will  be  expressed  in 
the  form  of  inequalities.  These  inequalities  will  be  the  first  conditions  that  we  need  to  assure 
stability. 

Step  5:  The  Be  formulated  Second  Variation 

The  second  variation  at  an  equilibrium  is  thus  of  the  form 

Z?’(H  +  C^)=  (square) 

—  N^'N  (f  r*  X  6m  da  +  f  St  x  m*  da) 

Jo  Jo 

'([  p*  X  6m da+  f  6r  x  m* da) 

Jo  Jo 

^rxm'ds)!!’ 

lla'H  Jo  Jo 

+  2^"(|la*||’)o*  ®  o*(  f  T*  X  Sxnda  +  f  6r  x  m*  da) 

Jo  Jo 

•  ( /  p*  x  6m  da+  I  Sr  x  m*  da) 

Jo  Jo 

—  2 - x-a  -ii  SrxSxnda)+l  — Sm'Smda+  I 

||tt«||>  Vo  '  Jo  Po  Jo  9a  da 

(3.16) 

Where  we  note  that 


N’'N  =  (-  ^^1+  2*»(||a*||»K  ««•) 


»  ra  I 


I  1^1  |  »»A  I 


fJ-1  _  “ 


I  +  2^«(lla‘l|*)tt*  ®  a*)-» (-  +  2^"(lla*ll»)a*  ®  a‘) 

a* 


=  QeJ.Q. 


(3.17) 


Step  8:  Collect  Integrals  of  Cross  Products 


Collecting  terms  containing  the  integrals  of  cross  products  the  second  variation  can  be 


written 


I>2(H  +  C*)=  {square) 

—  (QcJcQe  —  Qe)(  /  r*  X  £mds+  /  Srx  m‘ds) 

Jo  Jo 

■{I  r*  X  Sin  d8+  I  Sr  x  m*  ds) 

Jo  Jo 

ut‘  •  oi*  ,  f ^  1  dSr  dSr 

—  2 - r®*  *  (  /  ^rx5mds)+  /  — fm-fmd5+  /  K-^--^ds 

lla*ll*  Vo  Jo  Po  Jo  Be  ds 


(3.18) 


Step  7:  A  Vector  Identity 

Observe  that  a  simple  vector  identity  enables  us  to  write 

"  “■  ■  <r  = *  r*  w“' 

where  we  have  used  the  skew-symmetric  matrix  5(x)  associated  with  the  cross-product 

0  ~*S  *2 

S(x)  =  *s  0  -xi  ( 

— X2  Xj  0 

Step  8:  A  Quadratic  Form 


(3.19) 


(3.20) 


Now  define  the  symmetric  matrix 


R  =  Q.J.Q.-Q. 


(3.21) 


JS 


We  will  see  below,  that  an  eigenvalue  estimate  (S.23)  relies  on  having  R  nonnegative  definite. 
We  thus  require  that  conditions  on  the  parameters  of  the  problem  and  ^"(||a*||’)  bold  such 


that  defined  in  (3.13)  is  positive  definite  andR  defined  in  (3.21)  is  nonnegative  definite. 
The  latter  will  assure  that  R  has  a  square  root  We  will  examine  these  assumptions 


again  in  remark  2  below. 


Expanding  the  second  term  in  (3.18),  we  can  re-express  it  as  a  quadratic  form, 


R(/'  r*  X  Sm  da  +  f  Sr  x  m*  da)  •  {  f  r*  x  £m  da+  f  Sr  x  m*  da) 

Jo  Jo  Jo  Jo 

nR(S(r*(s))5in(s)  -  S(m*(a))5r(s))  •  (S(r*(a))^m(a)  -  S(m'(a))«r(cr))  deda 

_ 

=  _s^(m*(a))  -S(m*(a))]  dado 

^  jo  jo  «r^(«)]A^(a)A(a)  dado  (3.22) 


We  now  can  find  a  lower  bound  on  the  above.  The  bound  we  want  is  obtained  from  an 


eigenvalue  inequality  which  we  introduce  by  way  of  the  following  lemma 


Step  9:  An  Eigenvalue  Inequality 


Lemma(3.1):  Let  A(a)  €  L3^’*(0,/),  and  x(a)  €  LJ(0,^)  then 


nx^(a)A^(a)A(a)x(a)  dada  <  /  x^(a){  /  IA®(a)  da}x(a)  da,  (3.23) 
I  Jo  Jo 


where  A® (a)  is  the  maximum  eigenvalue  of  A®^(a)A(a). 

Proof:  Let  ||  •  ||  denote  the  standard  norm  in  Euclidean  space  and  also  the  induced  matrix 


norm  associated  with  it.  Then 


nx^(a)A^(a)A(a)x(<T)d<Tda<  f  f  |x^(a)A^(a)A(o)x(ff)| dads 
j  Jo  Jo 

^  f  f  l|A(»)x(a)||  |lA(<7)x(a)lldffda 
Jo  Jo 

<  r||A(a)||||x(a)||da  f' \\A{a)\\Mo)\\do 
Jo  Jo 


where  we  have  used  ||A(a)x(a)||  <  ||A(a)||  ||x(a)||.  We  can  now  use  the  Schwarz  inequality 


(/‘||A(.)||||x(.)||*)><  /‘||A(.)||»i.  /''||x(,)||»<, 
Jo  Jo  Jo 


Fin&lly  noting  that  the  value  of  ||A(«)||  is  simply  the  square  root  of  the  maximum  eigenvalue 
of  A^(s)A(s)  establishes  the  result.  | 

If  we  let  A^(s)  be  the  maximum  eigenvalue  of 

r  S^(r*(s))RS(r‘(s))  -S^(r«(s))RS(m«(s))-l 


A^(s)A(s)  = 


-S^(m*(«))RS(r*(s))  S^(m‘(«))RS(m'(s)) 


(3.24) 


and  let  =  /Q*A*(«)ds  then  we  have  by  way  of  lemma  3.1  a  lower  bound  on  the  second 


variation 


D*(H  +  C^)(p. >  {square) 

—  A*  /  6n 
Jo 


Sm^Sm  ds 


Jo 


■  /o  ®<  llo'IP 

+  f  ~6m^Smds+  f  (3.25) 

Jo  Po  Jo  08  ds 

Step  10:  A  Poincare  Type  Inequality 

If  we  assume  that  K  is  diagonal  and  use  a  Poincar4>type  inequality 

/^K^-^ds>e  ^K^r.^rds,  (3.26) 

Jo  V®  Jo 

with  c  =  then  the  second  variation  can  be  bounded  below  as 

£>’(H  +  C^)(p.,rMn«)  >  (  ) 

—  A*  /  Sroi^Sm  ds  —  A*  f  ds 

Jo  Jo 

~  Jo  ‘ 

+  [  Ism^Sindo  +  e  f  5r*'K5rrfs  (3-27) 

Jo  Po  Jo 

Step  11:  Rewrite  The  Lower  Bound 

We  can  reformulate  the  lower  bound  in  a  clearer  form  as  follows 

I?’(H +  C^)(p. >  {square)^ 

fi\  rtfml  r«mi 

+  /  -  fif  '  ^  (3.28) 

Jo  cK  -  IA»  J  L  J  L  J 


)  jr  ds 


If  we  define  the  matrix 


D(p%r%m*)  = 


ii-iP  -S(^»') 


(3.29) 


then  we  can  state  the  following  theorem  ; 

Theorezn(3.2):  If  the  matrix  R  =  QcJeQc  —  Qe  defined  in  (3.21)  exists  and  is  nonnegative 
definite,  J(  defined  in  (3.13)  is  positive  definite,  and  the  matrix  D  defined  in  equation  (3.29)  is 
positive  definite,  then  the  system  described  by  equations  (1.2)*  (1.4)  is  nonlinearly  (formally) 
stable  at  the  equilibrium  point  (p*,  r‘,  m*).  | 


Remark  1:  This  result  establishes  only  formal  stability,  since  it  is  based  on  the  the  definiteness 
of  second  variation.  To  establish  rigorous  stability  of  the  nonlinear  system  one  generally  needs 
to  examine  convexity  estimates  as  is  done  in  [2]. 

Remark  £:  Note  that  if  Q7*  exists  and  we  use  the  matrix  inversion  lemma  [4,  p.656]  we  obtnin 
the  following 


=  -R  (3.30) 

Recall  that  we  already  have  an  assumption  of  nonnegative  definiteness  on  R.  Thus  we  need 
to  specify  conditions  on  the  parameters  and  ^"(||a*||’)  such  that 

J"*  +  Qe  >  0  (3.31) 

(J  +  Qr*)"‘  <  0  (3-32) 


which  are  the  same  conditions  as  R  >  0  and  J«  >  0.  In  the  examples  of  the  next  section  Qe 
is  singular. 

Remark  S:  A  better  result  can  be  had  by  observing  that  A^(8)A(s)  is  frequently  in  the  form 
of  a  block  diagonal  matrix 


A*'(.)A(s)  = 


■Ar(.)A,(.) 

0 


0 

Ar(»)A*(«). 


(3.33) 


where  0  <  A:  <  6  and  because  of  the  semidefiniteness  of  A^A(s)  some  of  the  diagonal  blocks 
may  be  zero.  If  we  let  A?(s)  be  the  maximum  eigenvalue  of  A^(«)Ai(s),  0  <  t  <  A;  then  we 


can  define 


D'  = 


\ill 


(3.34) 


Thus,  if  the  conditions  of  theorem  (3.2)  are  satisfied  and  also  the  matrix  D'  defined  in  equation 
(3.34)  is  positive  definite,  then  the  system  described  by  equations  (l.2)-(l.4)  is  (formally) 
nonlinearly  stable  at  the  equilibrium  point  (p*,  r*,  m*).  In  theorem  (3.2)  this  will  mean  the 
special  choice  A*(s)  =  max{A}(s),...,A|(s)}. 

4.  Some  Examples 

In  this  section  we  apply  theorem  (3.2)  to  specific  equilibria  of  (l.l)  •  (1.3).  We  will  assume 
that  that  the  linear  extensible  shear  beam  lies  along  the  same  direction  as  the  second  principal 
axis  of  inertia  of  the  rigid  body.  From  geometric  considerations  the  position  of  the  shear  beam 
will  cause  the  principal  axes  of  the  rigid*body-shear>beam  configuration  to  lie  in  the  same 
directions  as  those  of  the  rigid  body.  In  this  case  the  addition  of  the  shear  beam  will  have 
the  efi'ect  of  increasing  the  moments  of  inertia  about  the  first  and  the  third  principal  axes. 
Because  the  linear  extensible  shear  beam  cannot  deflect  laterally  the  principal  axes  of  the  of 
the  configuration  remain  fixed  for  any  longitudinal  extension  of  the  shear  beam.  Thus,  for 
this  configuration  there  are  three  axes  about  which  the  equilibria  can  exist.  These  axes  will 
correspond  to  the  three  principal  axes  of  the  rigid  body. 

4.1.  A  Trivial  Equilibrium 

The  simplest  case  to  be  considered  is  when  the  rotation  takes  place  about  the  axis  along 
which  the  linear  extensible  shear  beam  lies.  In  this  case  the  equilibrium  will  be 

«•  =  W]  Sa  (4.1) 

r*  =  (aa  +  s)  €2,  0  <  s  <  €  (4.2) 

m*  =  0  (4.3) 


This  describes  the  linear-extensible-shear-beam  being  unstretched. 


What  follows  is  a  special  case  of  the  second  variation  computed  in  Step  1  of  the  previous 
section.  In  this  and  the  following  example  we  will  assume  ^"(||oc*||^)  is  the  same  as  in  [2],  thus 
recall  from  (3.7)  that  if  this  is  the  case  then 


2jMr 


And  the  two  quantities,  J,  and  Qe,  which  we  define  in  Step  4  are 


Q. 
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For  our  example,  if  we  first  compute 


=  122^2  «2 


then 

=  y2a(u;2)*  and  =  y|2(u;2)* 

from  which  we  immediately  compute 
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and  finally. 
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(4.4) 

(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


For  Jc  to  be  positive  definite  we  require  >22  >  Jiit  and  J22  >  Jss-  This  will  assure  positive 
elements  along  the  diagonal  in  the  inverse  above. 


15 


•J 


Thus,  the  quantity  Qe  JeQc  which  appears  in  the  reformulated  second  variation  of  Step  5 
be, 
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where  we  have  used  equation  (3.17)  and  the  fact  that  J*  and  Qe  are  diagonal. 

We  also  need  the  skew  symmetric  matrix  which  appears  in  Step  7.  Thus,  we  compute 

r  0  0  u;;1 

S(£f5"‘)  =  0  0  0  (4. 


l-w}  0  0 

Now  we  compute  R,  which  is  defined  in  Step  8. 


R  =  Q.J.q.  +  ^{l-^) 
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which,  along  with  the  definition  of  S(-)  in  (3.20),  we  can  now  use  to  compute 


S^(r*)RS(r*)  =  *  0*  0  0 

raa— Jii  * 


S^(r*)RS(in*)  =  0 
S^(m*)RS(m*)  =  0 


1%  ♦ 


These  matrices  are  used  to  form  the  matrix  A~(«)A(«)  in  (3.34),  note  that  it  has  only  the 
two  nonzero  elements  (computed  in  (4.15)).  These  correspond  to  the  first  and  second  diagonal 
elements.  Hence,  A^(«)A(s)  is  a  diagonal  matrix  and  the  nonzero  eigenvalues  are  these  two 
elements.  As  a  consequence  we  will  use  the  modified  bound  described  in  Remark  S.  Thus,  the 
eigenvalue  inequality  is  easily  obtained. 

After  using  the  Poincar^  inequality  of  Step  10 'nt  proceed  to  the  final  step  and  construct 
the  D'  matrix  in  (3.34) 


— - 1 — ^ —  Jn  rJ*  ds  0 


D'  = 


0  «- 5^77 /o' ’■5’* 


(a)'*. 


To  assure  that  the  D'  matrix  is  positive  definite  we  require 


J22  -  Jit  >  Po 


J22  -  J3S  >  Po 
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(4.17) 


(4.18) 


(4.19) 


and  also. 


(4.20) 


(4.21) 


Physically  the  first  two  conditions  are  classical  stability  conditions  on  the  stable  axes  of  rotation 
for  a  rigid  body.  The  term  on  the  right  is  the  additional  inertia  due  to  the  flexible  appendage 
which  adds  inertia  about  both  the  first  and  third  axes.  The  second  two  inequalities  are 
conditions  on  the  admissible  rotation  rates  of  the  configuration.  They  have  an  interesting 
physical  interpretation. 


4.2.  A  Non-TVivial  Equuibrium 


For  the  second  example  we  will  consider  rotations  of  the  rigid-body-shear-beam  configu¬ 
ration  about  the  first  or  third  principal  axes  of  inertia.  We  will  examine  the  case  when  the 
rotation  is  about  the  first  principal  axis  of  inertia,  rotations  about  the  third  axis  are  similar. 
This  corresponds  to  the  example  in  Krishnaprasad  and  Marsden  [2]. 
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(4.24) 


In  these  equations  we  have  0  <  s  <  f.  For  simplicity  we  will  denote  the  nonzero  element  of  r 
as  f],  and  that  of  m  as  m|. 

We  first  compute 

a*  =  +  f  r\m\d8kx  (4.25) 
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=  (inu;;  +  f  rlml  ds)wj 
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Subsequently  we  will  denote  the  first  element  of  a  by  ai.  We  now  compute 

a«o*^  fO  0  0‘ 

[o  0  ij 

and  finally,  and  Q*  defined  in  Step  4  ure 
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For  J,  ^  to  be  positive  definite  we  require 


tti  >  iiawi  oi  >  jssul 


(4.26) 


These  conditions  will  hold  if  jn  >  J321  <ind  jn  >  ja  and  will  assure  positive  elements  along 
the  diagonal  in  the  inverse  above.  These  conditions  are  the  same  as  (5.10)  in  [2]. 

Then  from  equation  (3.17)  we  have 
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The  skew  symmetric  matrix  of  Step  7  is 
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Now  we  compute  R  as  defined  in  Step  8, 


where  we  have 
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(4.31) 


(4.32) 
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(4.34) 


Note  that  these  are  not  the  same  as  the  71,  and  72  terms  which  appear  in  [2]. 


We  can  now  compute 


7ss»’2*  0  0 

S^(r')RS(r*)  =  0  0  0 

I  0  0  0 

S^(r‘)RS(m*)  =  0 
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From  this  we  can  compute  the  matrix  A^(s)A(s)  in  (3.24),  note  that  it  has  only  two  nonzero 
elements.  These  correspond  to  second  and  fourth  diagonal  elements.  Hence,  A^(s)A(s)  is 
a  diagonal  matrix  and  the  nonzero  eigenvalues  are  these  two  elements.  As  in  the  previous 
example  we  will  use  the  modified  bound  described  in  Remark  S. 

We  can  construct  the  D'  matrix  in  (3.34) 
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To  assure  that  the  D'  matrix  is  positive  definite  we  require 
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These  conditions  are  exactly  those  of  (5.14)  in  Krishnaprasad  and  Marsden  and  they  assure 
stability  about  the  equilibrium  which  also  satisfies  (4.26). 


Finally  a  remark  about  the  difference  between  [2]  and  our  development.  If  we  integrate  the 
matrix  we  call  A^(s)A(a)  then  the  elements  of  the  integrated  matrix  would  correspond  to  ‘72, 
and  '71  in  the  paper  of  Krishnaprasad  and  Marsden.  This  suggests  modifying  the  procedure  in 
the  previous  section  to  look  at  the  eigenvalues  of  the  integrated  matrix  rather  than  integrating 
the  eigenvalues. 
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